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ABSTRACT 



We report on a scheme for incorporating vertical radiative energy transport into a fully relativistic, Kerr-metric model of 
optically thick, advective, transonic alpha disks. Our code couples the radial and vertical equations of the accretion disk. 
The flux was computed in the diffusion approximation, and convection is included in the mixing-length approximation. 
We present the detailed structure of this "two-dimensional" slim-disk model for a = 0.01. We then calculated the 
emergent spectra integrated over the disk surface. The values of surface density, radial velocity, and the photospheric 
height for these models differ by 20%-30% from those obtained in the polytropic, height-averaged slim disk model 
considered previously. However, the emission profiles and the resulting spectra are quite similar for both types of 
models. The effective optical depth of the slim disk becomes lower than unity for high values of the alpha parameter 
and for high accretion rates. 

Key words, black holes physics — accretion disks 



1. Introduction 

Modeling accretion flows onto black holes is crucial for un- 
derstanding the energetic emissions observed from many 
sources, both Galactic and extragalactic (e.g., X-ray bina- 
ries, ULXs, AGNs). The accreting matter commonly settles 
into a disk-like configuration in which angular momentum 
is transported outwards as a result of shear in the differen- 
tially rotating fluid, allowing mass to be accreted onto the 
central compact object. Of particular interest are sources 
such as the low-mass X-ray binaries in which matter is 
transferred onto the X-ray source from a binary companion, 
usually a late-type star. When containing a black hole, these 
sources undergo irregular outbursts (lasting several months 
or years) that are analogous to the outbursts observed in 
dwarf novae, and are thought to be related to the viscous- 
thermal instability, which leads to an enhanced transfer of 
angular momentum in the accreting fluid (Lasota 2001). 

In outburst, the source is a luminous X-ray emitter 
and as the luminosity decays from the outburst maximum 
to minimum, the disk is thought to follow a sequence of 
quasi-stationary states, whose spectra can be derived from 
steady-state models of accretion disks. At high luminosities 
describing the inner regions of such accretion disks by the 
most-often used "thin disk" models is no longer valid. A 
description of such disks in terms of "slim disk" models is 
more appropriate. 

Following the approach of IShakura fc SunvaevI ()1973ll . 
in the "standard" discussion of stationary disk structure 



one assumes that there is no radial advection of heat (the 
flow is "radiatively efflcient"), the radial pressure gradi- 
ents are negligible, and the distribution of angular mo- 
mentum is Keplerian. All of these a ssumptions hold for 
very t hin disks, including the original IShakura fc SunvaevI 
(|l973f ) model (henceforth SS) and its general-relativistic 
ve rsion (NT) which was d eveloped for the Kerr geometry 
bv lNovikov Sz Thornd (|1973. 1. However, it has long been re- 
alized that not all disk-like accretion flo ws satisfy these as- 
sumptions In particular, slim disks fe.g.. lAbramowicz et al.l 
1988) form a stable branch of accretion, in which advection 
of entropy is important and the disk is not necessarily ge- 
ometrically thin. These solutions are particularly relevant 
to sources with high accretion rates, and they converge to 
the standard case in the low accretion rate limit. 

In the case of thin disks, methods for treating the ra- 
diative transfer and computing the vertical struc t ure o f 
the disk have been developed by Shaviv &: W ehrse (1986|) 
and|Hubenv. Q998i') and applie d by, e.g., Davis fc Hube;n vl 
(|2006D . IIdan et all (|2008D . and iRozanska fc Madeil (|2008l ). 



These have not yet been applied to slim disks. Until re- 
cently, the properties of slim disk models have been repre- 
sented with quantities averaged over the vertical structure 
of the disk. This paper presents improved advective steady- 
state models of accretion, in which the vertical structure of 
the "slim disk" is explicitly taken into account. 

Recently, ISq,dowskil ()2009D revisited sl im disk models 
with an improved numerical code, while ISq,dowski et al.l 
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l|2009 [) performed a preliminary study of the slim disk verti- 
cal structure. While our present work is based on these two 
papers, it breaks with the long tradition of computing the 
radial dependence of disk quantities, including the emergent 
flux F{r), before the vertical structure is analyzed. In this 
paper, by closely coupling the radial properties of the disk 
to its vertical structure, we offer a more consistent treat- 
m ent of the slim di s k. A similar approach has been taken 
bv lDotan fc Shavivl l|2010f) in their discussion of winds from 
Super-Eddington slim disks in a pseudo-Newtonian poten- 
tial. 

The high luminosities observed in the bright X-ray 
sources imply there is an effective mechanism of angular 
momentum transport in the accreting fluid, and much work 
has been devoted to MHD simulations of one such mech- 
anism, the magneto- rotational instability (MRI). For thin 
disks, the effectiveness of angular momentum transport is 
conveniently parametrized by the alpha parameter that was 
introduced by Shakura and Sunyaev (1973). Observations 
of systems in which variability is driven by accretion disk 
instabilities suggest rather high values of a for fully ion- 
ized disks, much higher than the value a = 0.001 that was 
adopted in the original slim-disk article (jAbramowicz et aD 
Ll988). Smak (1999) has already showed that one has to take 
a fti 0.2 to describ e hot disks in dwarf novae. As reviewed by 
iKing et all (|2007[ ). the value of a in hot, fully ionized disks 
has to be ~ 0.1 — 0.4, whereas in cold protostellar and FU 
Urionis disks, much lower values of this parameter are re- 
quired ~ 0.01 and ~ 0.001 — 0.003, respectively. Numerical 
simulations give a value of a ~ 0.01 for ionized disks, ten 
times lower than suggested by observations, and it proba- 
bly reflects the lim itations of the present MRI calculations 
([King et al.ll2007t) . However, the apparently observational 
determinations of a mentioned above are in fact strongly 
model-dependent. In the case of dwarf-novae, for example, 
these determinations assume that outbursts are described 
well by the thermal- viscous instability model. However, the 
case of the epitome dwarf-nova SS Cyg shows that such a 
description might be inadequate (jSchreiber fc Lasotall2007t 
ISmakll20Tof) . putting presumably observational determina- 
tions into doubt. Therefore, one cannot exclude that the 
numerical simulations are correct after all. 



In this paper we adopt the value a = 0.01. With such a 
low value of a, the resulting disk models are optically thick 
(tcS > 1), thus allowing a simplified treatment of radiative 
transfer. We treat the vertical energy transport in the dif- 
fusive approximation for the radiative flux (and account for 
vertical convection in the mixing- length approximation). 

We begin the body of the paper with a brief discussion of 
the model, then both the radial and the vertical equations 
of disk structure are presented. In Sect. [3] we describe the 
numerical method used to solve the problem. The radial 
structure of the solutions is described in Sect. 14. ll while the 
vertical profiles of various disk quantities are discussed in 
Sect. By comparing our results with previous work, we 
find that, with a judicious choice of parameters, a standard 
polytropic slim disk model can be made to approximate our 
full solution. Analytical formulae for the appropriate values 
of the parameters can be found in Sect. [SJ Finally, in the 
last section we discuss the limitations of the model and its 
possible applications. 



2. Slim-disk structure 

Slim disks are alpha disks. Like Shakura-Sunyaev disks they 
are based on the assumption that the dissipation mecha- 
nisms operating in accretion flows may be described by a 
viscous stress tensor whose leading component is propor- 
tional to the pressure. 

Slim disk solutions, again just like the SS and NT solu- 
tions, are obtained by solving vertically averaged (or height- 
integrated), radial equations of motion. Thus, steady slim 
disks, like the SS and NT ones, neglect the vertical struc- 
ture of flow, and describe essentially flat fluid configura- 
tions. Although an expression for the radial dependence 
of disk thickness is obtained, the slope of the disk surface 
is usually neglected in the discussion of emergent spectra, 
where a plane-parallel atmosphere is typically considered 
for the purposes of computing radiative transfer in the ver- 
tical direction. 

In the case of slim disks, the radial equations of struc- 
ture are a set of ordinary differential equations (ODEs). 
For a detailed account of these Kerr-metric equations 
(and state-of-th e-art traditional slim-disk models), see 
ISadowskil (I2009D. w ho followed previous work beginning 
with Lasotal (I1994D. and considered subsequent improve- 
ment s ( Abramowicz et aI1ll996l . Il997t ICammie fc PophamI 
Il998l ). The set of equations is closed by including rela- 
tions describing vertical hydrostatic equilibrium and ver- 
tical transport of energy. By solving these equations with 
appropriate boundary (or regularity) conditions, one ob- 
tains the radial profiles of the central temperature, Tc(r), 
surface density, S(r), the height-integrated pressure, P(r), 
the half-thickness of the disk, h{r), a radial velocity V{r), 
the emerging fiux of radiation, F{r), and certain other phys- 
ical quantities. 

In this paper we are taking a first step towards con- 
structing truly three-dimensional slim disk models, by solv- 
ing a set of differential equations describing the vertical 
structure of the disk. The resulting z-dependence of phys- 
ical quantities is used to compute certain coefficients that 
enter the radial equations and that up till now have been 
estimated with algebraic expressions. We refer to the result- 
ing models as "two-dimensional (2-D) slim disk" solutions, 
although it has to be understood that i n contrast to the 
thin-disk solutions of lUrpinl (|1984f ) and iKluzniak fc Kits] 
(2000); the actual meridional flow has not been computed 
in this paper, but the average radial velocity alone has been 
considered in the structure equations. However, the models 
are now self-consistent in the sense that the vertical aver- 
ages of physical quantities that form the coefficients of the 
radial ODEs do correspond to the vertical structure consid- 
ered in the radiative transfer calculation — in previous work 
the vertical structure of the radiative atmosphere was con- 
sidered a posteriori, and it had no influence on the radial 
structure of the disk. 



2.1. Basic assumptions, parameters, and coefficients 

We assume an axially symmetric, stationary fluid configu- 
ration in the Kerr metric, with fixed values of the black hole 
(BH) mass, M, and spin, a, parameters and the fluid disk 
is symmetric under reflection in the equatorial plane of the 
metric. Matter is supplied at a steady rate, M, through a 
boundary "at inflnity" and angular momentum is removed 
through the same boundary (in practice, we use the NT 
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solution for the outer boundary condition), whereas zero 
torque is assumed at the BH horizon. We assume that no 
mass or angular momentum crosses the disk surface. We 
neglect the loss of angular momentum to both wind and 
radiation. We do not consider self-irradiation of the disk 
and assume that the magnetic pressure may be neglected. 
Neglecting the incoming radiation may not be justified for 
super-Eddington accretion rates for which the disk is geo- 
metrically thick. 

A fraction. (1 -I- f^'^^{r))^^, of the entropy generated 
locally by dissipative processes is released into the radiation 
field, while the remainder is advected by the gas. 

A unique solution to the slim-disk model can only be 
found if certain additional assumptions are made. We make 
the following arbitrary choice. We neglect the vertical vari- 
ation {z dependence) of the velocity field, considering only 
its height-averaged value; thus, the velocity is always di- 
rected radially inwards and is a function of the radial co- 
ordinate alone. Similarly, we assume there is no z variation 
of the advection factor f'^'^^{r). Dissipation and angular 
momentum transport are given by the alpha prescription 
([Shakura fc Sunvaevlll973[ ). with a constant value of a. We 
assume that the dissipation rate is proportional to the total 
pressure, p. For a more detailed statement, see Eq. (|13p and 
the comment following it. Calculations are carried out for 
the value a = 0.01. 

We are looking for 2-D slim-disk solutions at a definite 
value of mass accretion rate for a given Kerr metric. Thus, 
for a fixed value of a, there are three fundamental param- 
eters describing a given slim-disk solution: A/, a, and M . 

In the structure equations, we take G — c = 1 and make 
use of the following expressions involving the BH spin: 



A 
A 
C 
V 
H 



- 2Mr + a, 
-^r^(? + 2Mra'^, 

0„ ^-3/2 



= 1 - 3r; 
= l-2r: 



with a, — a/M and = r/M. 

We use the Boyer-Lindquist system of coordinates and 
introduce the vertical coordinate z = r cos9. The metric 
near the equatorial plane takes the form (jAbramowicz et aD 

[Mi, 



— —df + —(d(j)- ujdtf + —dr^ 



dz^ 



(3) 



where w = 2Mar/A is the angular velocity of the frame 
dragging. 

The radial gas (three-)velocity, V, as measured by an 
observer co-rotating with the fluid at a fixed v alue of r, is 
given by the relation (jAbramowicz et aLlllQOGf ) 



v/Vi - 



ru 

Ai72^ 



(4) 



where is the contravariant radial component of the 
fluid 4-velocity and has the dimension of physical ve- 
locity. The disk surface density is E = J^j^ pdz, while 

the height-integrated pressure is P = J^j^pdz. The to- 
tal pressure is the sum of gas and radiation pressures. 



P = Pgas +Prad- We adopt an equation of state correspond- 
ing to the choice pgas = kpT j {^.m-p), and prad = aT*/3, 
with k the Boltzmann constant, rrip the proton mass, and 
a the radiation constant (no confusion with the spin pa- 
rameter may arise) . The mean molecular weight is taken to 
be /i = 0.62, but see the comment following Eq. (HI]). The 
height-integrated energy density is 



E 



h 

-h 



Ps 



1- 



1 



3pr 



ad 



(5) 



where 7 = 5/3. 

The following averages (moments) enter the radial equa- 
tions as coefflcients of certain terms: 



m 

V3 



+ h 



1 

JO 

2 

^ Tc Jo 
E/P, 

i r 



dz, 

-h 

pT dz, 



pz"^ dz. 



(6) 

(7) 
(8) 
(9) 



All integrals in this section are taken at a flxed value of r. 
Here, T{z) is the gas temperature, and Tc is its value at the 
equatorial plane, r(0) = Tc- 

The vertical epicyclic frequency squared, which c an be 
thoug ht of as the vertical component of gravity, is (jKatol 
11995 



(1) We also define 



02 _ MH 



ri = i + — . 



(10) 

(11) 



(2) 2.2. Vertical structure equations 



We describe the vertical structure of an accretion disk in 
the optically thick regime by th e following equat ions: 
(i) Hydrostatic equilibrium (jKato et al.ll2008l ). 



p dz 



(12) 



where n± is defined in Eq. (|TU)) . Although other expressions 
for the right- hand side of Eq. (| 1 2|) can be found in the lit- 
erature (e.g., lAbramowicz "et aLni997t ) includes 7^ 0), 
the form above is appropriate for our scheme, in which 
the vertical structure is precalculated before any informa- 
tion about the radial variables becomes available. Thus, in 
Eqs. (|10p and (jl2p . we assume Keplerian angular velocity 
[M/(Cr3)]i/2 _ as wen as Vz = in hydrostatic equilibrium. 

(ii) The energy generation equation. We assume that 
the vertical fiux of energy inside the disk is generated 
according to 



dz 



32? 
2C 



ap 



1 + / 



1/2 



(13) 



Strictly speaking, this does not correspond to a constant 
a prescription, as the term {5/2){M/r^y^^ , which is de- 
rived from Keplerian strain, departs somewhat from the 
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value that would follow from the actually computed distri- 
bution of angular momentum (cf., Fig. [ij. However, as the 
departure for sub-Eddingtonian accretion rates is small, we 
expect Eq. (|13p to afford a good approximation. Note that 
jadv _ Q corresponds to the NT disk (going over into the 
Shakura-Sunyaev disk in the nonrelativistic limit of thin 
disks), /^'^^ > 1 characterizes advection-dominated disks, 
while /^'^^ < describes those disk regions where the ad- 
vected heat is being released. The amount of heat advected 

(iii) Energy transport. The structure of the disk has to 
be such that the actual value of the divergence of the flux 
corresponds to Eq. (jl3p . Radiative transport is computed 
in the diffusive approximation 



J-(z) 



16crT3 dT 
Znp dz ' 



(14) 



while convective transport is computed in a mixing-length 
approximation. Energy is transported in the vertical direc- 
tion through diffusion of radiation or convection according 
to the value of the thermodynamical gradient, which can 
be either radiative or convective. Accordingly, we take 



dlnT 
dlnp 



Vrad, for Vrad < Vad 
V conv ; 

for Viad > Vad 



(15) 



The thermodynamical quantities Cp, Vad and 
(dhi pI d \n.T)r, are calculate d using standard formu- 
lae (e.g.. IChandrasekhaij|1967[ ) assuming solar abundances 
{X = 0.70, Y = 0.28) and, when necessary, taking the 
effect of partial ionization of gas on the gas mean molecular 
weight into account . 

We use Rosseland mean opacities (including 
the processes o f abso rptio n and scatte r ing) taken from 
[Alexander et all (I1983D an d'Seato n et all (|1994[ ). Following 
other authors (e.g., Idan e t al. 200^, we neglect expansion 
opacities, in agreement with our neglect of vertical velocity 
gradients. 

(iv) We set the following boundary conditions. At the 
equatorial plane (z — 0) we set — 0, in accordance 

with the assumption of reflection symmetry, while at the 
disk surface, T{h) = 0, we follow the Eddington approxima- 
tion (|Mihalaslll98ll and require F{h) = aT^g = 2(jT^{h). 

In practice, for a fixed r, and prescribed values of Tc and 
Z^"^^, a trial value of the central density, pc, is assumed and 
the equations are integrated in z until /^(z,) = 10~^^ gcm~'^ 
(as a stand-in for the disk surface, = h). If J-{h) and T{h) 
fail to satisfy the surface boundary condition, the assumed 
value of pc is adjusted, and the integration is repeated, 
until the condition J^{h) = 2crT^(/i) is met. Convergence is 
usually attained in a few iterations. The emergent flux of 
radiation at any given r is then F — T{h). 



with the adiabatic gradient given by a derivative at con- 
stant entropy: Vad = {d\nT/dhip)s- 

The radiative gradient Vrad is calculated in the diffusive 
approximation. 



v.. 



fad 



3pK rF 
16crT4rj2 



(16) 



where k/; is the Rosseland mean opacity, and a is the 
Stefan-Boltzmann constant. At the equatorial plane we ap- 
ply the boundary conditions described in point (iv) below. 

When the temperature gradient exceeds the value of 
the adiabatic gradient, we have to consider the convective 
energy flux. The convective gradient Vconv is calculated 
using the mixing length theory introduced by iPaczvhskil 
l|1969D . We take the following mixing length. 



Hrn\ — 1.0i/„ 



(17) 



with pressure scale height Hp defined as (jHameurv et al.l 
[19981) 

= pnlz + ^pn^- ^^^^ 

The convective gradient is defined by the formula 

^conv — ^ad + (V,.ad-Vad)j/(j/ + «;) (19) 

where y is the solution of the equation 
9 tI 



y + wy + w y — w = 0, 



(20) 



with the typical optical depth for convection Tmi 
pKRHmh and w given by 



3r„ 



bl2a^T<^Hp V^lnT 



(V,,d - Vad) ■ (21) 



2.3. Radial structure equations 

The radial sector of the model is described by four laws of 
conservation and a regularity condition: 
(i) Mass conservation. 



M ^ -27rSru''. 
(ii) Conservation of angular momentum, 

M ^ , ^i/2Ai/2r 

— (£-An) = aP, 

Ztt r 



(22) 



(23) 



where C = is the specific angular momentum, is a 
constant, whos e value is to be specified la ter, and F is the 
Lorentz factor (jGammie fc PoDhamlll998h : 



F2 = 



1 



£2^2 



i-y2 ' A ' 

(iii) Conservation of radial momentum, 

V dV _A ]_dP 

1 — dr r Y. dr ' 

where 

MA {n-n+){n -Q-) 



1 - n^R^ 



(24) 



(25) 



and n = u'^/u'^ is the angular velocity with respect to a sta- 
tionary observer, f2 = fi — is the angular velocity with re- 
spect to an inertial observer, ilf = ±M^/2/(j.3/2 -|-ajVfi/2) 
are the angular frequencies of the corotating and counter- 
rotating Keplerian orbits and R = A/{r'^A^^^) is the radius 
of gyration. 

The value of P is taken from vertical structure solutions 
for given values of E and Tc , 



P 



prrip 



(26) 
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hence, the radial derivative of the height-integrated pres- 
sure takes the form 

dr dr dr 



dr dr 



(27) 



with /3 = ri2{k / ptm-p)YiTc/ P ■ 
(iv) Energy conservation. 
The advective cooling is defined following iKato et "aTl 
([200I) in terms of the vertically integrated quantities, as 

Id c\P rin 

Q^''" ^-i-iru^{E + P))-u^—- u'-^dz. (28) 
r dr dr J oz 

Using mass conservation (Eq. [22]) and hydrostatic equi- 
librium (Eq. [12]), the expression can be rewritten as 



^adv 



M ( P d In P 



'73 



27rr2 V S dlnr 



(l + '?3) 



PdlnS 
E dlnr 



P d In 773 
E dlnr 



dln?74 
dlnr 



(29) 



Just like P, the advective cooling term, Q'*'*^, and the coef- 
ficients ?73, ?74 are all determined from the vertical structure 
solutions. 

The manipulation of the last term in Eq. (pS)) was as 
follows: 



dp 

- / u'^—dz^n] u'^pzdz^ 
J-h dz J_,, 

27rr dr 

3. NUMERICAL METHOD 

3.1. Vertical structure 

The set of ordinary differential equations describing the ver- 
tical structure, i.e., Eqs. p^ . and ([T5]) together with 
appropriate boundary conditions are solved for a given BH 
spin on a three-dimensional grid spanned by the radius r, 
the central temperature Tc, and the advection factor f^'^'^. 
For a given set of these parameters, we start the integration 
from the equatorial plane (z = 0), and the solution satis- 
fying the outer boundary condition is found as described 
at the end of t)2.2l The resulting quantities describing the 
vertical structure (Tc, S, Q^*^^, P, 771, 772, 773, and 774), to- 
gether with r, are printed out to tables for subsequent use 
in interpolation routines. As it turns out, the first two of 
these parameters, Tc and E, can be used to uniquely de- 
termine all the other quantities characterizing the vertical 
structure, including f^*^^ (see Fig.[2|). 

Calculating the full grid of vertical solutions for a single 
value of BH spin takes about 5 hours on a 4-CPU worksta- 
tion. 



3.2. Radial structure 

By a series of algebraic manipulations of Eqs. ([Hj), 

, , and , we obtain the following set of two or- 
dinary differential equations for V{r) and Tc{r) 



1 dlnV 
1 - y2 dlnr 



-Do' 



(30) 



(4-3/3) 



dlnPc 
dlnr 
27rr2 E 



■Da 



ladv 



r(r - M) 
A 

din 771 



(/3-ri) 

^dln772 



d In r d In r 
d In 773 E 774 d In 774 



dlnr 

with A/i and Vq given by 



P 773 d In r 



(31) 



Ml ^ A 



Mm 

din 773 
dlnr 



Q 



adv 



p 
s 



r(r - M) 



2 774 d In 774 
773 d In r 



-Do 



(32) 



(33) 



Typically, Vq vanishes close to the black hole, as C ap- 
proaches dn- 

To obtain a solution one has to solve this system of two 
ordinary differential equations, together with the following 
regularity conditions at the sonic radius rs , defined by the 
same conditions; 



Do\ 



0, 



(34) 



as well as outer boundary conditions given at some large ra- 
dius rout ■ The solution between the outer boundary a nd the 
sonic point is found using the relaxation technique (jPressI 
l2002f ). with Cin treated as the eigenvalue of the problem. 
The sonic point for a* = is located at 5.9M for accre- 
tion rate 0.01 AlEdd and at 5.0M for 2.0MEdd- To start the 
relaxation process we have to provide a trial solution that 
is obtained by a method simil ar to the one descr i bed in 
ISadowskil (pOOOl) assuming the iNovikov k Thornd ifTOT^ 
outer boundary conditions. Once the trial solution is found, 
one can start the relaxation process with a free inner bound- 
ary corresponding to the location of the sonic point. To find 
the solution inside the sonic point we make a small step in- 
ward to cross the critical point and then integrate down 
to BH horizon using a Runge-Kutta method of the fourth 
order. In this work we use 25 mesh points spaced logarith- 
mically in the radius on the section between the sonic point 
and rout = lOOOAf . This particular number of grid points 
is enough to resolve all disk features. We have verified that 
the results are accurately reproduced with a denser grid. 

The parameters linked to the vertical structure (P, 
Q'^'^", rji, 772, 773, and 774) for given E and Tc are hn- 
early interpolated from pre-calculated tables of the ver- 
tical structure solutions (for any value of V, S is deter- 
mined directly from mass conservation, Eqs. (j22p l. The 
radial derivatives din 771 /dlnr, dln 772/dlnr, dln 773/dlnr, 
and dln 774/dlnr are evaluated numerically from the 7^1, 772, 
773, and 774 profiles in the previous iteration step. A relaxed 
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4 6 8 10 12 14 16 18 20 



r [M] 

Fig. 1. Angular momentum (u^) in a Schwarzschild slim 
disk for two accretion rates. For very low accretion rates the 
angular momentum follows the Keplerian profile (dotted 
line) down to the ISCO. For high accretion rates the flow 
is super-Keplerian between the "center" of the disk at rcon 
and the "potential spout" at rpot- The vertical dot-dashed 
line on this and subsequent figures denotes the location of 
the ISCO. 

solution is obtained in a few iteration steps. Once a solution 
outside the sonic point is found we numerically estimate the 
radial derivatives of V and Tc at the sonic point using val- 
ues given at r > rs and use these derivatives to start direct 
integration inside the sonic point. 

The solution thus obtained may then be used as a trial 
solution when looking for the relaxation solution of another 
slim disk, i.e., when one of the three fundamental param- 
eters (M, a, M) has a slightly different value. Each relax- 
ation step takes approximately 5 seconds on a single-CPU 
workstation. 

4. RESULTS 

In the following two sections we present and discuss both 
the radial and vertical structure of slim accretion disks. 
All the solutions, if not stated otherwise, were computed 
assuming a = 0.01 and M = 10 Mq. 

4.1. Disk radial structure 
Angular momentum 

The angular momentum profiles for our accretion disk so- 
lutions near a nonrotating BH are presented in Fig. [T] 
Results for two values of mass accretion rate are shown, 

M = O.lMEdd and 2.0AfEdd, where AfEdd = 16LEdd/c^ is 
the critical accretion rate that for a disk around a non- 
rotating BH approximately corresponds to the Eddington 
luminosity, -Z^Edd- For the lowest accretion rates the profiles 
follow the Keplerian profile and reach its minimal value 
(£in in Eg. (1^) a t the marginally stable orbit (ISCO, 
iBardeen et al., ,19721 ). The higher the accretion rate, the 
stronger the deviation from the Keplerian profile. The disk 
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Fig. 2. The T^-S plane at r = 20M for a nonrotating BH 
(a* = 0). The dotted lines connect solutions for the ver- 
tical structure of slim disks that have the same value of 
the advection parameter f^'^^. The locus of standard (ra- 
diatively efficient, f^'^^ = 0) disk solutions is shown with 
the thick dotted line. The solid thick line represents the 
vertical slim-disk solutions for different accretion rates (in- 
dicated by triangles), and the dashed line presents corre- 
sponding solutions of the conventional polytropic slim-disk 
model (see The difference between the two lines in the 
low M limit corresponds to the difference in E between the 
two models (see Fig. [TTll . 



is sub-Keplerian at large distances and super-Keplerian at 
moderate radii. The Keplerian profile is crossed again at a 
point located inside the marginally stable orbit, and corre- 
sponding to what is usually called "the cusp" or "the poten- 
tial spout". For a detailed study of the physics of the inner 
edge of a see lAbramowicz et aU (|2010( ) . 

S-curves 

Figure [2] presents slim disk solutions at r = 20M on the 
Tc-S plane, for a nonrotating BH. Solutions of the poly- 
tropic, height-averaged models are presented for compari- 
son; for detailed discussion see ^ The locus of solutions 
for various values of the m ass accretion rate ha s the shape 
of the so-called "S-curve" (EE ramowicz et al.llTQSH'l . The 
lower, gas-pressure dominated branch accurately follows 
the track of radiatively efficient solutions (/'^'^^ = 0). The 
middle, radiation-pressure dominated branch is reached at 
M w O.lMEdd- As advection becomes significant, the slim- 
disk solution leaves the f^'^'^ = track and moves to higher 
advection rates (S-curve). Around M = 5MEdd, the solu- 
tions enter the upper advection-dominated branch corre- 
sponding to /^'^^ > 1.0 (more than 50% of heat stored in 
the accreted gas). At M = 20A/Edd this rate almost in- 
creases up to 80% iP'^'' w 4.0). 

Surface density 

Profiles of the surface density for 2-D slim-disk solutions 
for a non-rotating BH are presented in the upper panel 
of Fig. [3) Different regimes, corresponding to different 
branches of the "S-curve" on the (S, Tc) plane are visi- 
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Fig. 3. Profiles of the surface density (upper panel) and 
corresponding values of radial velocity V (bottom panel) 
of a slim disk for a nonrotating BH. Solutions for different 
accretion rates are presented. 



ble. For large radii the surface density increases with in- 
creasing accretion rate (the lower gas-pressure dominated 
branch), while this relation is opposite for moderate radii 
(the middle radiation-pressure dominated branch) . For ac- 
cretion rates M > 5.0MEdd the upper advection-dominated 
branch would be reached. The local maxi ma in the surface 
density profiles (discussed in detail in, e.g.. lSq,dowskil[2009h 
are visible for moderate accretion rates (^ O.SMEdd)- The 
bottom panel of Fig. [3] presents corresponding profiles of 
the radial velocity V as measured by an observer corotat- 
ing with the fluid. 

The surface density dependence on BH rotation is pre- 
sented in Fig.m The profiles are shifted to lower radii as the 
inner edge of the disk moves inward for higher BH spins. 
The outer parts of the accretion disk are insensitive to the 
metric. 



Optical depth 

In the top panel of Fig. [S]we plot the total optical depth of 
the vertical slim disk solutions for different accretion rates 
(a* — 0). The total optical depth 



10" 



1000 




Fig. 4. Profiles of the surface density for slim disks at a 
constant accretion rate [M = O.lMEdd) and various BH 
spins. 



where the total opacity coefficient kh, which includes the 
processes of absorption and scattering, is closely related to 
the surface density. Indeed, the radial profiles of the optical 
depth shown in Fig. [5] follow the corresponding profiles of 
surface density. Any differences in the profiles come from 
the dependence of the opacity coefficient on local density 
and temperature. Outside the ISCO the total optical depth 
is always large {rtot > 10^). 

The diffusive approximation for radiative transport may 
only be used if photons are absorbed, otherwise LTE can- 
not be established. In a scattering-dominated atmosphere, 
the effective optical depth is then the relevant quantity to 
be used in checking for the self-consistency of the diffusive 
approximation. The bottom panel of Fig. [S] presents cor- 
responding profiles of the effective optical depth, which is 
estimated in the following way: 



Vi 



,)K[i pdz 



Ttot 



krp dz. 



where k^s — 0.34 cm'^g ^ is the mean opacity for Thomson 
electron scattering. For M > 0.3A/Edd the inner region 
of the disk becomes radiation-pressure dominated (it en- 
ters the middle branch on the corresponding S-curve), the 
surface density decreases with increasing accretion rate, 
and electron scattering begins to dominate absorption. 
Therefore, the effective optical depth decreases with in- 
creasing accretion rate and reaches values < 1 in the 
inner parts of the disk for accretion rates above l.OMEdd- 
As a result, for M > l.OA/Edd, the diffusive approximation 
can no longer be applied, and our present approach to solv- 
ing for the disk structure breaks down. In Fig. [6] we exhibit 
the dependence of the effective optical depth on the value 
of a. In general, t^s is inversely proportional to a (as is the 
surface density). At lower accretion rates (M < O.lMEdd) 
the effective optical depth remains large even for high val- 
ues of a. 
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Fig. 5. Optical depth for a — 0.01 slim disks around a 
Schwarzschild black hole at different accretion rates. Top 
panel: The total optical depth as a function of the radius. 
Bottom panel: The effective optical depth as a function of 
the radius. The ISCO is shown at r = 6M. The shaded 
region in the bottom plot indicates the region where the 
diffusive approximation is invalid. 
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Fig. 6. Profiles of the effective optical depth of a 
Schwarzschild slim disk for three values of viscosity (a — 
0.01, 0.05, and 0.1), calculated for two accretion rates, 
O.lMEdd (dashed lines) and l.OAfEdd (solid lines). 
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Fig. 7. Flux emitted from the surface of a slim disk at five 
accretion rates onto a Schwarzschild black hole. At high 
accretion rates significant emmision from within the ISCO 
is clearly visible in the figure. 



Flux profiles 



Profiles of the flux emitted from the disk surface {F = 
cT^ff) in the case of a nonrotating BH are presented in 
Fig. [T] Results corresponding to accretion rates from 0.01 
up to 5.0MEdd are shown. For the lowest rates the emis- 
sion from inside the marginally stable orbit is negligible 
as expected in the standard accretion disk models. This 
is no longer true for higher accretion rates, and the ad- 
vectio n of energy causes signifi cant emission from smaller 
radn (jAbramowicz et all [20To[) . For super-Eddington ac- 
cretion rates the emitted flux continues to grow with a 
decreasing radius even inside the marginally stable orbit. 
Radiation coming from the direct vicinity of the black hole 
is suppressed by the gravitational redshift (the (/-factor). 
Therefore, an observer at inflnity will observe a maximum 
in the profile of the effective temperature even for the high- 
est accretion rates. 

The increase in the advective flux with increasing accre- 
tion rate is clearly visible in Fig. [51 The ratio of the heat 
advected to the amount of energy emitted is presented for 
different accretion rates. For very low accretion rates these 
profiles approach the limit of a radiatively efficient disk 
^yadv ^ Q-j^ advection component becomes signif- 

icant for higher accretion rates. Some part (up to 30% at 
r = 20Af for 2.0AfEdd) of the energy generated at mod- 
erate radii is advected with matter and radiated away at 
r < lOA/. This causes the significant change in the emitted 
flux profile at the higher accretion rates visible in Fig. [T] 

In Fig. [5] we present the emitted flux profiles for dif- 
ferent BH angular momenta at a constant accretion rate 
M = O.IA/Edd- These profiles coincide at large radii where 
the influence of the BH rotation is negligible; however, the 
higher the BH spin, the closer to the horizon the marginally 
stable orbit. Therefore, in the case of rotating BHs, the ac- 
creting matter can move much deeper into the gravity well, 
compared to nonrotating BHs. This effect leads to an in- 
crease in the disk luminosity and hardening of its spectrum. 
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Fig. 8. Profiles of the advection coefficient /^'^^ for differ- Fig. 10. The height of the photosphere at different accre- 
ent accretion rates (Schwarzschild black hole). tion rates onto a Schwarzschild black hole. 
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Fig. 9. Flux profiles at fixed accretion rate (M = O.IA/Edd) ^^S- Profiles of photospheric height at a constant ac- 



for five values of BH spin. 



cretion rate (M = O.lAfEdd) and various BH spins. 



which can be inferred from Fig. ^ — the higher the spin, the 
higher the disk luminosity, and the higher the flux (which 
corresponds to the effective temperature). 

Photosphere location 

The flux observed at infinity may be obtained by perform- 
ing ray tracing of photons emitted from the accretion disk 
(see ^ . Scattering in the layer above the photosphere must 
also be taken into account. An accurate calculation requires 
detailed knowledge of atmospheric properties, including the 
location of the photosphere. This is particularly important 
when accretion rat es are high and the disk is no longer 
geometrically thin (jSe|,dowski et al.ir2009[ ). 

In Fig. [TO] we plot the profiles for the z-location of the 
photosphere, i/phot, obtained in our model at different ac- 



cretion rates for a nonrotating black hole. Clearly, for high 
accretion rates, M > O.lMEdd, the inner regions become 
thicker (effects of radiation pressure). For M = I.OAfEdd, 
the ratio ffphot/f reaches a value as high as 0.25. Close to 
the ISCO the height of the photosphere rapidly decreases 
because of vigorous cooling (compare Fig. Although 
the rapid change in disk thickness violates the assumption 
of the hydrostatic equilibrium that we make when solv- 
ing for the disk vertical structure, the accelerations con- 
nected with the vertical motions involved are much lower 
than the vertical component of gravity. Indeed, the vertical 
accelerations are close to Vrdvz/dr ~ VrdivrdH / dr) / dr ~ 
vld{H/r)/dr - r^\{H/rf - ^\H{H/r). In deriving this 
estimate we liberally assumed that dH/dr ~ 1 and used the 
fact that the rapid decrease in disk height occurs near the 
sonic point, while the speed of sound Vs is approximately 
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Fig. 12. Comparison of the dynamical {vrdvz/dr « 
Vd{VdH/dr)/dr) and gravitational {n\H) components of 
the vertical equilibrium equation at the photosphere. The 
stars denote locations of the sonic radii. 



r^l± (H/r) . Thus the acceleration terms modifying Eq. 
would be smaller than the gravitational acceleration by a 
factor of a few percent: {H/r) ^ 10~^. In Fig. [T^ we plot 
the dynamical and gravitational components of the verti- 
cal equilibrium equation at the photosphere. It is clear that 
the former is at least 10 times smaller than the latter at the 
sonic radii for M < MEdd- Therefore, in all likelihood, our 
results correctly describe the disk structure that would be 
obtained without assuming strict hydrostatic equilibrium. 

In Fig. [TT] we present radial photosphere profiles, at a 
fixed accretion rate and different values of the BH spin. 
The photospheric heights coincide for large radii. In the 
inner regions of the disk, the height of the photosphere 
increases with BH spin, reflecting the increased luminosity 
and radiation pressure. 

4.2. Vertical structure 

In Fig. [13] we present a vertical cross-section of a 
Schwarzschild slim disk for M = O.OlMEdd- At this accre- 
tion rate the disk is radiatively efficient and no advection 
of entropy is expected. The top panel presents the radial 
profiles of the photosphere and the disk surface (defined as 
a layer with p = 10~^^g/cm^). 

The total optical depth reaches values as high as ~ 
20000 on the equatorial plane and decreases monotonically 
towards the photosphere at r = 2/3 (Fig. [121 second panel 
from the top). Within the marginally stable orbit, the total 
optical depth significantly decreases, as shown in Fig. [S] 

The third panel of Fig. [T^] presents the temperature, 
which decreases with height from at the equatorial plane 
to Teff at the photosphere. The maximum value of the tem- 
perature is attained on the equatorial plane at r ~ lOM. 

Density is presented in the fourth panel, and the next 
panel presents the vertical flux that is generated inside the 
disk according to Eq. [T31 It is set to zero on the equatorial 
plane by reflection symmetry, and then rapidly increases, 
because in an alpha disk the dissipation is proportional to 




Fig. 13. Vertical structure of a slim disk for M = O.OlAfEdd 
and = 0. The top panel presents the surface of the disk 
(green dashed line) , and the photospheric surface (red solid 
line). The other panels present the structure of the disk 
below the photosphere. Top to bottom: total optical thick- 
ness, temperature, density, vertical flux of energy, and the 
termodynamical gradient. The blue solid line in the bottom 
panel delimits the convective region. 



the pressure, which reaches its maximum on the z — 
plane. Close to the disk surface, where pressure is almost 
negligible, the flux slowly settles down to the emitted value. 
At the accretion rate chosen for the flgure the flux rapidly 
decreases inside the marginally stable orbit. 

The bottom panel of Fig, [represents the nonmonotonic 
distribution of the termodynamical gradient (Eq. ([TS]) '). 
which ranges between 0.2 and 0.4. Therefore, the disk's ver- 
tical structure cannot be described by a simple polytropic 
relation. Moreover, in the region of the highest temperature 
(close to the equatorial plane at moderate radii, r w 20Af ), 
the heat is transported upward through convection. 

The vertical structure of an accretion disk with ten 
times higher accretion rate, M = 0.1 MEdd, is presented in 
Fig- [HI The general picture remains the same, since the ad- 
vection of heat is still insignificant. However, as the temper- 
atures increase, the inner disk regions become dominated 
by radiation pressure. For O.lMEdd the convective region 
extends from the marginally stable orbit up to 300Af and 
covers more than half of the disk thickness. 
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The disk structure is significantly different in the case 
of a high accretion rate (e.g., l.OA/Edd), with a significant 
amount of advection. The vertical cross-sections of the slim 
disk are presented in Fig. 1151 The inner regions are domi- 
nated by radiation pressure, so the disk geometrically thick- 
ens and the photosphere is now higher. The total optical 
depth mostly follows the surface density (compare Fig. [SJ , 
and therefore decreases considerably towards the black hole 
in the inner parts of the disk. 

The temperature maximum is again located at the equa- 
torial plane close to r ss lOM. The maximum of the effec- 
tive temperature (corresponding to the vertical flux shown 
in the fifth panel (compare also Fig. [7]) is shifted inwards, 
down to r « 8M. 

The fourth panel of Fig. [13] presents the density distri- 
bution. Despite the fact that the surface density monoton- 
ically increases outwards (see Fig. [3]) , p has two maxima in 
the equatorial plane: at r 200Af and r 6M. Finally, the 
bottom panel presents the termodynamical gradient distri- 
bution. For such a high accretion rate, the convective zone 
extends nearly to the photosphere for r < 200A/ and is 
present up to r = 2000Af . 

In Fig.[TH]we present the density in the meridional plane 
for three different accretion rates. The equatorial plane lies 
in the middle of each plot. The violet boundaries denote the 
photosphere. The two maxima of density at the equatorial 



plane for M — l.OA/Edd) are clearly visible also in this 
representation, as are some other features discussed above. 



5. Comparison with height-averaged slim disk 
solutions 

In this section we compare our 2-D slim disk solutions, 
where the radial structure equations are coupled to those 
for the vertical structure, with the standard polytropic slim 
disk model, in which the slim-disk equations an d properties 
are a veraged over the thickness of the disk fe.g.. lKato et al.l 
l2008l ). We note here that the slim disk solutions presented 
in one of our previous papers (|Sadowskl,2009) did not fol- 
low the polytropic formalism, assuming different relations 
between the vertically integrated quantities and their val- 
ues on the equatorial plane (e.g., S = '2-poH instead of 
Y, — 2InPoH^ where /^v is defined in Eq. ([37])) In this work, 
we also use a more general form of the energy equa tion 
(compare our Eq. jH]) with Eq. (6) in lSadowski (j2009f )V 

Let us assume the polytropic equation of state with the 
polytropic index N: p — Kp^^^^^ . The vertical integration 
of the hydrostatic equilibrium formula, Eq. (jl2p . gives 

/ 2 \ ^ 



11 



A. S^dowski et al.: Vertical structure of slim disks 




0.03 
OOl 




Fig. 16. Meridional profiles of density for three accretion 
rates: M = 0.01 (top), 0.1 (middle) and l.OMEdd (bottom 
panel) in (r,z) coordinates. The violet boundaries show the 
location of the photosphere. The black hole is described by 
Mbh = 10 Mq, and a* = 0. 



We also assume 



T = Tn 1 - 



(36) 



One can now calculate analytical formulae for 771 to 774 
(Sect. 



^0 "'0 



H 



/ TUz = hH 



H 



—— / pT dz 



Vi = 
where 



1 kl 



N+l 



J:Tc + 2haT^H 



N 



pz dz — JnH 



In 



J. 



r(l + iV) NeN {2^Nl) 
2 r(3/2 + A^) 

ir(3/2 + 7V) 



{2N - 
1 



4r(5/2 + 7V) 6 + 4iV 



(37) 



(38) 



In this approach we do not solve the vertical structure 
consistently, so we need to make some additional assump- 
tions about the vertical equilibrium of forces and radiation 
transfer. Following other authors, we simplify the hydro- 
static equilibrium (Eq. ((T^)) by applying a finite difference 
approximation and write 



2r)2 



(2iV + 3)|. 



(39) 



One has to remember that disk thickness, H, defined in this 
way is not the exact location of the photosphere. Therefore 
we introduce a factor fn relating these quantities. 



phot 



(40) 



Assuming that radiation is transported in the vertical di- 
rection through diffusion, the radiative flux is given by 



16crT3 dT 
3kp dz 



(41) 



Under the one-zone approximation one obtains the follow- 
ing formula for the total flux emitted from disk surface. 



(42) 



where factor fp has been introduced to account for inac- 
curacies arising from this approximation, as well as from 
the dominance of the disk convection in certain regions (as 
discussed in H4.2p . Now, advective cooling takes the form 



r3 



3SK 



(43) 



Usually, the following values of N, Jh and fp are as- 
sumed: 



N = 3.0 
Ih - 1.0 
fp = 1.0 



(44) 



In Fig. [T7] we compare radial profiles of the flux, pho- 
tospheric height and surface density of the 2-D slim disk 
model described in this paper (consistently taking its ver- 
tical structure into account) with profiles obtained for the 
polytropic, height-averaged, slim disk, as presented in this 
section with parameters defined in Eq. (HH). The com- 
parison was carried out for two accretion rates (0.1 and 
l.OMEdd), with a — 0.01. For the lower accretion rate 
(O.lMEdd), the disk is radiatively efficient (advective cool- 
ing is negligible) , and therefore both flux profiles almost co- 
incide and correspond to the Novikov & Thorne solutions. 
However, the disk photosphere location in the polytropic 
model turns out to be overestimated by more than 20% in 
the region corresponding to the maximal emission. The pro- 
files of the surface density do not coincide either, the 2-D 
solutions giving values ~ 25% lower than the corresponding 
polytropic solutions (compare Fig. [5]) . 

For the higher accretion rate (l.OMEdd), the flux proflles 
remain similar (up to 1%). However, as advection becomes 
important, the emission is shifted inwards with respect to 
the Novikov & Thorne profile. The photosphere location in 
the polytropic model is overestimated by ~ 30% and the 
surface density by ~ 20%. 

A question arises as to whether such differences in the 
fiux, photosphere, and surface density profiles affect the 
resulting disk spectrum. In Fig.[T5]we present spectral pro- 
files and their ratios (2-D to polytropic) for two accretion 
rates. The spectra were calculated with ray- tracing rou- 
tines Hursa 2006) using the BHSPEC package (iDavis et al.l 
l2005h . assuming the inclination angle i = 70° and dis- 
tance to the observer d = lOkpc. As BHSPEC gives tab- 
ulated solutions of the full, frequency-dependent, radiative 
transfer equations for the disk vertical structure taking the 
Compton scatterings into account in the disk atmosphere, 
the spectra presented in Fig. [TS] are not those of a simple 
multi-color blackbody. However, one should be aware that 
using BHSPEC for calculating spectral color correction is 
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Fig. 17. Comparison of the flux, disk thickness (the height 
of the photosphere for models presented in this paper and 
of the zero-density surface for 1-D polytropic models) and 
surface density profiles calculated using the 2-D (this paper, 
solid lines) and the usual polytropic {N = 3, dotted lines) 
slim disk models for a = 0.1. The solutions for two accretion 
rates (0.1 and l.OAfEdd) are presented with thick and thin 
lines, respectively. 

not consistent with the assumed vertical structure (calcu- 
lated or height-averaged) because it is based on a stand- 
alone disk model. 

The general shape of the spectra is similar for both types 
of slim disk models, because the emission profiles nearly co- 
incide. However, the spectra are not identical. The bottom 
panel of Fig. [TH] presents ratios of the spectral profiles of 
the corresponding solutions for each accretion rate. They 
all coincide att low energies (< O.lkeV), while for higher 
energies the discrepancies are as large as 3% for l.OAfEdd 
at 5keV. These differences are atributed to (slightly) dif- 
ferent profiles of the flux, the photosphere, and the surface 
density. We conclude that the proper treatment of the ver- 
tical structure hardly affects the spectra of slim disks for 
this range of accretion rates (M < MEdd) and the viscosity 
parameter (a < 0.01). 




energy [keV] 

Fig. 18. The upper panel presents spectral profiles of the 
2-D and polytropic solutions for two accretion rates (0.1 
and l.OA/Edd), at inclination angle i = 70° and distance 
d = lOkpc. The bottom panel presents ratios of the corre- 
sponding spectra (of 2-D to polytropic solutions) for both 
accretion rates. 

We end this section by giving fitting formulae for iV, Jh, 
and fp, which approximate the full numerical 2-D solutions 
described in this paper with a polytropic slim disk model. 
An advantage of using these formulae lies in avoiding the 
need to perform time-consuming calculations of the vertical 
structure and avoiding numerical problems connected to 
interpolation in the vertical solutions grid. The formulae 
for the polytropic model parameters for a = 0.01 are 

N = 3.25 X Sn 

f„ - 0.63 X Sh (45) 
Jf = 0.94 X Sf, 

where the spin correction coefficients Sn, Sh, and Sf are 
given by 

Sn = 1 + 0.002 {6-rrns/M), 
Sh = 1 + 0.003 (6 - r„,/M), 
Sf = 1 + 0.064 {6~r,ns/M). 

Here, rms is the radius of the marginally stable orbit. For 
a nonrotating BH, at the radius r = 7M (corresponding to 
the highest disk effective temperature), the fitting formulae 
are accurate to 1% for the emitted flux, the photospheric 
height and the surface density. 

6. Discussion 

Motivated by a desire to explain and fit the observed spec- 
tra of accreting black holes in binary systems to theo- 
retical models, we have developed a 2-D model of opti- 
cally thick slim disks. These should be particularly rele- 
vant to transient binaries. In quiescence, their inner disk 
regions are described well by optically thi n advection- 
dominated accretion flows (ADAFs; see e.g., iLasota et al.l 
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llQQGtlDubus et aT1l200lf ). However, a few black hole systems 
(e.g., GRS 1915+105 and LMC X-3) have been observed 
in thermal sta tes corresponding to disk luminosities higher 
than O.SLEdd (|McClintock fc Remillardll2003l ; ISteiner et all 
l2010t) . and modeling these require optically thick models 
going beyond the standard thin disks. 

In this work we present a 2-D slim disk model, in which 
the radial and vertical structures are coupled. Such an ap- 
proach eliminates arbitrary factors that influence solutions 
of the usual polytropic slim disk model. The results were 
obtained under two key assumptions: an alpha disk was 
assumed (dissipation proportional to pressure) , with a uni- 
form value of a, and the fraction of the generated entropy 
that is advected was computed at every radius under the 
assumption that this fraction does not vary with the height 
above the disk plane (Eq. [T3]). Both of these assumptions 
seem arbitrary, and we can offer no physical motivation for 
the (conventional) choice we made. 

Under these assumptions and for the value a — 0.01 
of the viscosity parameter, we computed and presented 
the detailed structure of 2-D slim disks, parametrized by 
the mass accretion rate, and the two Kerr metric param- 
eters, M and a. Somewhat surprisingly, the spectra ob- 
served at infinity from such disks differ by only a few per- 
cent from those obtained from previously considered slim 
disk models (in which the equations and structure corre- 
spond to a height average over a polytropic atmosphere). 
Such differences are unlikely to introduce any large correc- 
tions to spin measurements based on X-ray continuum fits 
made with corresponding height-averaged polytropic mod- 
els of slim disks. However, already the latter produce sig- 
nificantly softer spectra in the sub-Eddington regime than 
the Novikov & Thorne model. For high luminosities, fits 
based on slim disk models may therefore provide higher 
values of the black hole spin parameter than corresponding 
fits based on the kerrbb model (e.g., Shafee et al. 200(|). 
This issue will be discussed i n detail in a forthcoming pa- 
per (jBursa fc Sadowskill20H"l ). 

One has to be aware that the model of vertical structure 
presented here is only an approximation of the real phys- 
ical processes taking place in disk interiors. The diffusion 
approximation and the convection treatment in the mixing 
length approach are known to successfully describe media 
with large effective optical depths but break down when the 
disk becomes optically thin. We have shown that the effec- 
tive optical depth of slim accretion disks may drop below 
unity for super-Eddington luminosities and sufficiently high 
values of a. For such conditions, a more sophisticated model 
of radiation transfer should be implemented. However, for 
a < 0.01 and L < L^dd the assumptions of this work are 
self-consistent. For higher values of a their range of appli- 
cability is limited to lower luminosities (e.g., to 0.5-LEdd for 
a = 0.1). Kerr-metric slim disks with low eff ective opti- 
cal depth were discussed by iBeloborodovl (|l998f ) , who finds 
them to be significantly hotter than the optically thick ones. 
We have already started implementing a radiative transfer 
scheme valid for disks with small effective optical depths 
into the scheme introduced in this work. It will be presented 
and discussed in a future paper. 

Another remark is connected to the fact that one can 
expect winds to be blown out of the disk surface at super- 
Eddington luminosities. Such a phenomenon may signifi- 
cantly change the disk structure, e.g., its thickness. This 



feature of slim disks, not describe d in our calculations, ha s 
recently been cleverly modeled bv lDotan &: ShavTvl (j2010[ ). 
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